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Abstract: 

We continue our study of the construction of analytical coefficients of the epsilon-expansion 
of hypergeometric functions and their connection with Feynman diagrams. In this paper, 
we apply the approach of obtaining iterated solutions to the differential equations 
associated with hypergeometric functions to prove the following result: 

Theorem 1: 

The epsilon-expansion of a generalized hypergeometric function with integer values of pa- 
rameters, 

pFp-i{Ii -\- aie, ■ • • ,Ip + apE] /p+i -|- hiS, • ■ ■ , hp-i + z) , 

is expressible in terms of generalized polylogarithms with coefficients that are ratios of 
polynomials. 

The method used in this proof provides an efficient algorithm for calculating of the higher- 
order coefficients of Laurent expansion. 
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1. Introduction 

Hypergeometric functions are useful in the evaluation of Feynman diagrams. See, for ex- 
ample, Ref. H] for a review of how these functions arise. In this paper, we will be concerned 
with the manipulation of hypergeometric functions ^, by which we understand specifi- 
cally 

(1) the reduction of the original function to a minimal set of basis functions, 

(2) the construction of the all-order e-expansion of basis functions. 

The e-expansion refers to the Laurent expansion of hypergeometric functions about rational 
values of their parameters in terms of known functions or perhaps new types of functions. 
In the latter case, the problem remains to identify the full set of functions which must be 
invented to construct this expansion for general values of the parameters.^ 

Problem (1) is a purely mathematical one. It is closely related with the existence 
of algebraic relations between a few hypergeometric functions with values of parameters 
differing by an integer, the so-called "contiguous relations" The systematic proce- 
dure for solving the relevant recursion relation is based on the Grobner basis technique. 
In particular, a proper solution for generalized hypergeometric functions, the so-called 
"differential reduction algorithm," was developed by Takayama[^. (See Ref. for a 
review.) By a differential reduction algorithm, we will understand a relation of the type 
F{a lb j, b; z) = n'^-|^D(a + k; b)F(a, b; z), where j, k are integers, 6 is a list of additional 

^AU these procedures coincide with standard techniques used in the analytical calculation of Feynman 
diagrams. ^, ^ It has long been expected that all Feynman diagrams can be represented by some class of 
hypergeometric functions. Now we can propose specifically that any Feynman diagrams can be associated 
with the Gelfand-Karpanov-Zelevinskii (GKZ or A- hypergeometric function) hypergeometric functions 
Let us recall that Lauricella's, Horns' and generalized hypergeometric functions occur as special cases of 
the GKZ-systems. For an introduction, we recommended Ref. M. 
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parameters and D is a differential operator of the form D = A{a; b;z)-j^ + B{a; b; z)? For 
Gauss hypergeometric functions, the reduction algorithm was presented in Refs. |11, 12|. 
For generalized hypergeometric function, is it equivalent to the statement that any function 
pFp-i{a;b; z) can be expressed as a hnear combination of functions with arguments that 
differ from the original ones by an integer, pFp_i{a + fh; b + k; z) , and the function's first 
p — 1 derivatives: 

Rp+i{a,b,z) pFp^i{a + m;b + k; z) = 

i i?i(a, 6, z) ^-^^ H ^ Rp-i{a,b,z) ^ + Rp{a,b,z)\pFp_i{a;b;z) , (1.1) 



where m, k are hsts of integers, the Ri are polynomials in parameters at, bj and z. 

Problem (2) arises in physics in the context of the analytical calculation of Feynman 
diagrams. The complete solution of this problem is still open. We will mention here some 
results in this direction derived by physicists. Let us recall that there are three different 
ways to describe special functions: 

(i) as an integral of the Euler or Mellin-Barnes type, 

(ii) by a series whose coefficients satisfy certain recurrence relations, 

(iii) as a solution of a system of differential and difference equations (holonomic 
approach) . 

For functions of a single variable, all of these representations are equivalent, but some prop- 
erties of the function may be more evident in one representation than another. These three 
different representations have led physicists to three different approaches to developing the 
e-expansion of hypergeometric functions. 

The Euler integral representation (i) was developed intensively by Davydychev, Tarasov 
and their collaborators ||l^ and the most impressive result was the construction of the all- 
order e-expansion of Gauss hypergeometric functions in terms of Nielsen polylogarithms 



1 14]. This type of Gauss hypergeometric function is related to one- loop propagator- type 
diagrams with arbitrary masses and momenta, two-loop bubble diagrams with arbitrary 
masses, and one-loop massless vertex-type diagrams. 

The series representation (ii) is a very popular and intensively studied approach. The 



first results of this type were derived by David Broadhurst [15| for the so-called "single 
scale" diagrams, which are associated with on-shell calculations in QED or QCD, with 
further developments appearing subsequently in Ref. |l^.'^ Particularly impressive results 
involving series representations were derived recently by Moch, Uwer, and Weinzierl in the 



^An algorithm of differential reduction of generalized hypergeometric functions to a minimal set allows 
the calculation of any Feynman diagram that is expressible in terms of hypergeometric functions without 
any reference to integration by parts or the differential equation technique. The application of this algorithm 
to the calculation of Feynman diagrams will be presented in another publication. 

^Some relations between the Mellin-Barnes representations and series representations of Feynman dia- 
grams follow from the Smirnov-Tausk approach 
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framework of the nested sum approach. [p!7|, 15] Algorithms based on this approach have 



been implemented in computer code.[l^] The series approach leads to algebraic relations 
between the analytic coefficients of the e-expansion, but does not provide a way to obtain a 
reduction of the original hyper geometric function before expansion. Another limitation of 
this approach is that the parameters are restricted to integer values or special combinations 
of half-integer values (so-called "zero-balance" parameter sets).^ 

For approach (iii), obtaining iterated solutions to the proper differential equations 
associated with hypergeometric functions, the first results were obtained for Gauss hyper- 



geometric functions expanded about integer values of parameters. [ 22 1 In Ref. |^], that 
result was extended to combinations of integer and half-integer values of parameters. An 
advantage of the iterated solution approach over the series approach is that it provides 
a more efficient way to calculate each order of the e-expansion, since it relates each new 
term to the previously derived terms rather than having to work with an increasingly large 
collection of independent sums at each subsequent order. 

The aim of the present paper is to apply approach (iii) to proving the following theo- 
rem:^ 

Theorem 1 

The all-order e-expansion of a generalized hypergeometric function pFp-i{A + ae; B + be; z), 
where A and B are lists of integers, are expressible in terms of generalized poly logarithms 
(see Eq. \l.^ ) with coefficients that are ratios of polynomials. 

To be specific, this means that 

P{{a},{b},z) pFp-i{A + ae;B + be;z) = "^Rg{z) Ug{z)e'' , 

where s = (si,--- ,si) is a multiple index and P {{a} , {b} , z) , Rg{z) are polynomials. The 
generalized polylogarithms are defined by the equation 

Lifci,fc2,...,fc„ (z) = 2_. fci k2 fcT ' ^^-2) 



For completeness, we recall that generalized polylogarithms ( |1.2| ) can be expressed as iter- 
ated integrals of the form 

dt dt dt dt dt dt dt dt . . 

Lii., ... I. (z) = / — o O • • • O O 0---0 O 0---0 — o , (1.3) 

^' • " ^ ^ in i i t 1-t t t t 1 - 1 ' ^ ^ 

fci-1 times k„-i times 

where, by definition 

dt dt dt dt _ ["^ dti [^^ dt2 [^'^'^ dtk^_i dtj,^ 



o — o • • • o 

s.'t t ^ 1 — t Jq tl Jq t2 Jq tfei-l Jq 1 — ifci 

times 



*For some new results on e-expansions of hypergeometric functions with nonzero-balance parameter sets 
of parameters (specifically, one half-integer parameter), see Ref. po[ ]. 

^In fact, the result expressed in this theorem can be proved within the nested sum approach How- 
ever, our idea is to extend the iterated solution approach to this more complicated system and in the 
process, derive a more efficient algorithm for calculating the analytical coefficients of the £— expansion. 



-3- 



The integral (|1.3| ) is an iterated Chen integral |27] w.r.t. the differential forms ujq = dz/z 
and oJi = so that 



iJn 



UJl - ■ ■ Wq" UJl 



(1.5) 



2. All-order e-expansion of generalized hypergeometric functions with in- 
teger values of parameters 

In this section, we shall prove Theorem 1. We begin by noting that Eq. ( |1.1[ ) can be 
written in a slightly different form: in terms of any basic function pFp-i{a; b; z) and its 
first p — 1 derivatives, 

Rp+i {a, b, z) {a + rh;b + k] z) = 

|i?i(a, 6, z) QP-^ + ■■■ + Rp-i{a, b,z)9 + Rp{a, b, z)|pFp_i(a; b; z) , (2.1) 

where m,k are lists of integers, the Ri are polynomials in parameters ai,bj and z, and 
9 = z-^ . The essential step in proving Theorem 1 is the following lemma: 

Lemma 1 

The all-order e-expansion of the function pFp-.i(ae;l + be;z), is expressible in terms of 
generalized polylogarithms (Eq. 



Lemma 1 could be proved in the same manner as in case of multiple (inverse) binomial 
sums. This was done in Ref. ||l^. However, it is fruitful to prove it using the construction 
of an iterated solution of the proper differential equation related to the hypergeometric 
function.^ We will follow this technique here, and in the process construct an iterative 
algorithm determining the analytical coefficients of the epsilon expansion. 

Let us consider the differential equation for the hypergeometric function uj{z) = 
pFp_i(ae; 1 + be; z): 



p-i 



:Y[{e + a,e)-el[{9 + b,e) 



4 = 1 



1=1 



uj{z) = 



(2.2) 



The boundary conditions for basis functions are w(0) = 1 and 9^uj{z)\_^_^ = 0, where 
j = 1, • • • ,p — 1. The proper differential equation for uj{z) is valid in each order of e. 
Defining the coefficients functions Wk{z) at each order by 



(2.3) 



fc=0 



^The proper solution for Gauss hypergeometric functions was constructed in Refs. pi] , psf . 
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the boundary conditions for the coefficient functions are 



wo{z) 

Wk{z) 

{z^J-'wkiz) 



2=0 



2=0 



A; < , 

k>l , 

k>0 , 

k>0 . 



(2.4a) 
(2.4b) 
(2.4c) 
(2.4d) 

(2.4e) 



The differential equation (2^) has the form 



(l-z) 



dz 



dz 



p-i 



p-i 



Wk{z) = Yj 



i=l 



Pi{d)--Qi{h) 



dz 



p—i 



Wk-i{z) + Pp{d)wk-p{z) , 



where Pj{a) and Qj{b) are polynomials of order j defined on spaces of p- and (p— l)-vectors 
a and b, respectively. They are defined as 



Po = Qo = 1 , 

Pr = ... ,;^^]^nj^<...<j^aj^ • • • Oj, 



r = 1 



5 ;f ) 



Hi, ^...^iM, ■■■bi^ , r = 1, • • • ,p - 1 , 



Qr — ^ii,... ,ir=l'^-^ii<---<ir"«i ■ ■ ■ ) 

Qp = 0, 



(2.5a) 
(2.5b) 
(2.5c) 
(2.5d) 



so that 



Pi - Sj=iaj , 
P2 = ^Ij^i-i^jaiaj , 

-^3 = ^ri,i2,i3 = l;n<i2<i3"H^«2'^«3 ' 



P -U^ a- 



Qi = ^z[b, , 



Q2 = ^j^i-i^jbib 



3 ' 



Q3 = sr^ 



n,«2,«3=l;n<«2<«3^*i^*2^j3 



(2.6) 



Qp-i = ^^i=lbi , 
Qp = 0. 



The polynomials Pj and Qj satisfy the following relations: 

Pj{a, b) = Pj{d) + 6P,_i(a) , Qj{d, b) = Qj{d) + 6Q,-i(a), j = 1, • • • ,P • (2.7) 
In particular, 

Pj (a, 0) = Pj (a) , Qj (a, 0) = Qj (a) . 
Let us introduce a set of a new functions p^^\z),j = 1, • • • ,p — 1 defined by 



pO){z)=e^u{z)^[z-) io{z) = Y,PX\^)e\ j = l,---,p-l, 



(2.8) 



fc=0 
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where the coefficient functions p\^\z) satisfy 



The boundary conditions for these new functions follow from Eq. ( 2.4e| ): 

/5?(0) = 0, k>0, j>l. 
Eq. ( |2.5| ) can be rewritten as a system of first-order differential equations 

z^/^\z) = p^^^'\z) , j = 0,l,...,p-l 



i=l 



and we have 



Pi{a)-^Qi{b) 



Wk{z) = Pk\z) . 



The solution of system can be presented in an iterated form: 



dt 
1-t 



ftp it) 



4 = 1 

^Qi{b)p'iZi~^\z) - Qp_i{b)[wk-p+i{z) -6o^k-p+i] 



i=l 



jp'iKt), k>l, i = l,2, 



,p-l 



where 6a^b is the Kronecker delta function. 

From the system of Eq. ( |2.13D , it is easy to find that 

p'j^\z)=0, k<p; j = 0,l,--- 



(2.10) 

(2.11a) 
(2.11b) 

(2.12) 



(2.13a) 
(2.13b) 



(2.14) 



The first nonzero terms are generated by Eq. (2.13a) for i = k = p. Substituting this result 
into Eq. ( 2.13b ) we will find the solution of the first iteration: 

>-i-j)(z) = Pp{a)Ui+j (z) , j = 0, 1, • • • ,P - 1, (2.15) 



where Lij (z) is a classical polylogarithm |24| and Lii (z) = — ln(l — z). Lemma 1 fol- 
lows from the representation ( ^.llb ), the value wo{z) = 1, the definition of generalized 
polylogarithms ([L^), and Eq. ( |2.15D . 

The case when one of the upper parameters of the hypergeometric function is a positive 
integer number, pFp^i^Ii, A + ae; B + be;z), corresponds to ai equal to zero. A smooth 
limit exists in this case and the particular result can be reproduced from expression ( 2.13| ) . 
Theorem 1 is thus proved. 
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3. Explicit expressions for the first five coefficients of the expansion 



Let us return to Eq. ( ^.13| ) and look at the next terms of the expansion. The first nonzero 
term of the iteration is given by Eq. ( p.l4| ). The second iteration corresponds to k = p + 1. 
In the r.h.s. of Eq. ( p. 13a ), only terms with i = 1 produce a non-zero contribution. 



Aiiln2(l 
2 ^ 



QiLi2 (z) 



where for simplicity, we omit arguments in the functions Pj , Qj and introduce a notation 

Aj = Pj - Qj , j = 1, • • • , p - 1. 
Substituting the results in Eq. ( 2.13b| ), we will get the solution of the second iteration: 



AiLij+1,1 {z)-QiU2+j {z) , j = 0, 1, ■ ■ ■ ,p - 1, 



(3.1) 



where Lia^_... (z) is a generalized poly logarithm. The third iteration corresponds to 
k = p + 2, and in the r.h.s. of Eq. (2.13a) only terms with i = 1,2 will produce a non-zero 
contribution, 



AfLij+1,1,1 (z) + (A2 - QiAi)Lij+i,2 {i 



+ (Q?-Q2)Lij+3(z)-QiAiLij+2,i(z) , 
The fourth iteration corresponds to A; = p -|- 3 and equal to 



0,l,---,p-l. (3.2) 



P. 



AfLij+1,1,1,1 (z) + Ai (A2 - QiAi) [Lij-+i,i,2 (z) + Lij+1,2,1 (z)] 

+ (AiQ2 _ A1Q2 - A2Q1 + A3) Lij+1,3 (z) - QiA2Lij+2,i,i (z) 
+Qi (AiQi - A2) Lij+2,2 (z) + Ai {Qj - Q2) Lij+3,1 (z) 
- {Ql - 2Q1Q2 + Q3) Li,+4 (z) , J = 0, 1, • • • ,p - 1. 
The fifth iteration corresponds to k = p + 4: and equal to 



(3.3) 



P^^-t'\^) 

Pr, 



AfLij+1,1,1,1,1 (z) + {AiQi - 2A1A2Q1 + A|) Lij+1,2,2 (z) 

+AI (A2 - QiAi) [Lij+1,1,1,2 {z) + Lij+i, 1,2,1 {z) + Lij+1,2,1,1 (z)] 
+Ai {Ai {Ql - Q2) - A2Q1 + A3} [Lij+1,1,3 (z) + Lij+1,3,1 (z)] 
-QiA^Lij+2,1,1,1 (z) + QiAi (AiQi - A2) [Lij+2,1,2 (z) + Uj+2,2,1 (z)] 
+Al {Ql - Q2) Lij+3,1,1 (z) - Ai {Ql - 2Q1Q2 + Qs) Lij+4,1 (z) 

+Ql {Ai (Q2 - Ql) + A2Q1 - A3} Li,+2,3 {z) 

+ [Ql {Ai (Q2 - Ql) + A2Q1} - Q2A2] Lij+3,2 {z) 
+ {Qf - 3QIQ2 + 2Q1Q3 + Ql- Qa) Lij+5 {z) 

+ { A4 - Q1A3 + A2 {Ql - Q2) - Ai {Ql - 2Q1Q2 + Q3) } Lij+1,4 {z) , 
j = 0,1,--- ,p-l. (3.4) 
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For lower values of the index p, the following relations can be used for transforming har- 
monic poly logarithms to the classical |24] or Nielsen |25] ones: 



^,1,1, 



,1 



Sj-i,p+i ( 



p times 



So,,(^) = ^^ln^(l-^) 



Lii,2 (z) 
Lii,3 (z) 
Lii,4 (z) 
Li2,2 (z) 

Li3,2 (z) 

Li2,3 (z) 



ln(l-z)Li2(z)-2Si,2(^) , 
ln(l-z)Li3 (z)-^[Li2 (z)]^ 
ln(l - z)Lu {z)+F2{z) , 



^[Li2(z)]'-2S2,2(Z) , 

hi2{z)Lh{z) + ^F2{z)-2Ss,2{z) , 
1 



F2(z)--Li2(z)Li3(z) . 



1 



Lii,i,2 (z) = -ln^(l - z)Li2 (z) + 21n(l - z)Si,2(z) + 3Si,3(z) 
Lii,2,i (^;) = - ln(l - z)Si,2(z) - 381,3(2;) , 
Li2,2,i (z) + Li2,i,2 (z) = Fi{z) - U2 (z) 81,2(2;) , 



1 



1 



Lii,i,3 (z) + Lii,3,i (z) = ^ ln2(l - z)Li3 (z) + ^ ln(l - z) [U2 (z) 
- ln(l - z)82,2 (z) - Li2 (z) 81,2 (z) + Fi (z) 

Lii,2,2 (^) = ln(l - z) |2S2,2(^) - I [Li2 (z)]' 

+2Li2 (z) 8i,2(z)-2Fi(z) 

1 



Lii,i,i,2 (^) + Lii,i,2,i (z) + Lii,2,i,i (z) 



^ln3(l-z)Li2 (z) 





-- ln\l - z)8i,2(z) - ln(l - z)8i,3(z) - 2Si,4(z) 



(3.5a) 

(3.5b) 
(3.5c) 
(3.5d) 
(3.5e) 
(3.5f) 

(3.5g) 

(3.5h) 

(3.5i) 
(3.5j) 
(3.5k) 

(3.51) 
(3.5m) 
(3.5n) 



where we have introduced two new functions related algebraically (see Eqs. (2.23) - (2.25) 
in Ref. H): 



Fiiz) 



dx 



ln2(l -x)Li2 (x) 



fir 

F2{z) = / — ln(l-x)Li3(x) . 

For completeness, we will present the values of P and Q for p = 3, 4: 
• p = 2> 

Pi (a) = ai + 02 + 03 , Qi{h) = 61 + 62 , 

P2{a) = aia2 + 0103 + 0203 , Q2(p) = ■ 
P3{a) = 010203 , Qsib) = . 



(3.6) 
(3.7) 



(3.8) 
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• p = A 



Pi{a) = 01+02+03+04 , Qi{b) = 61+62+63 , 

P2{a) = 0102+0103+0104+0203 + 0204 + 0304 , Q2{b) = 6162 + 6163+6263 . 

^3(0) = 01O2O3+01O2O4 + 02O3O4 , Qsib) = 616263 . 

^4(0) = 01O2O3O4 , <54(6) = . (3.9) 



The first few coefficients, up to order 4, could be cross-checked using the results of Ref. 

Hi- 



We would like to point out that Eqs. ( p.lj ) - (3.4) contain an explicit logarithmic singu- 
larity at 2; = 1. It is well-known that the generalized hypergeometric function pFp_i(a; 6; z) 
converges absolutely on the unit circle \z\ = 1 if 



Re 



In this case, the coefficients of the e-expansion also converge at each order in e. To get a 
smooth limit, it is enough to rewrite Eqs. ( |3.1[) - ( |3.4[ ) in terms of functions of argument 
1 — z and set z = 1. 



4. Conclusions 



We have shown (Theorem 1) that the e-expansions of generalized hypergeometric func- 
tions with integer values of parameters are expressible in terms of generalized polyloga- 
rithms (see Eq. (|l.2|)) with coefficients that are ratios of polynomials. The proof includes 
(i) the differential reduction algorithm; and (ii) iterative algorithms for calculating the an- 
alytical coefficients of the e-expansion of basic hypergeometric functions (see Eq. ( |2.13| )). 
The first five coefficients of the e-expansion for basis hypergeometric functions are calcu- 
lated explicitly in Eqs. {^J^, (|]|), (pj), and (|3j). The FORM |9| representations 
of these expressions and the next coefficients are available via Ref. |30|. 
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